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ABSTRACT 

The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed 
cosmic star formation rate. This allows to estimate the rate expected in the horizon of the 
gravitational wave detectors advanced Virgo and ad LIGO and to compare those rates with 
independent predictions. In this context, the rapid neutron-capture process, or r-process, can 
be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic 
process. We compute the early cosmic evolution of a typical r-process element. Europium. 
Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach 
allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated 
evolution of Eu. 

We find that the bulk of Eu observations at [Ee/H] > -2.5 can be rather well fitted by 
either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution 
favors NSM as the main astrophysical site for the r-process. A comparison between our cal¬ 
culations and spectroscopic observations at very low metallicities allows us to constrain the 
coalescence timescale in the NSM scenario to ~ 0.1-0.2 Gyr. These values are in agreement 
with the coalescence timescales of some observed binary pulsars. Einally, the cosmic evolu¬ 
tion of Eu is used to put constraints on i) the NSM rate, ii) the merger rate in the horizon of 
the gravitational wave detectors advanced Virgo/ad LIGO, as well as iii) the expected rate of 
electromagnetic counterparts to mergers (’’kilonovae”) in large near-infrared surveys. 

Key words: Cosmology: dark ages, first stars. Physical Data and Processes: gravitational 
waves, nucleosynthesis. Stars: supernovae, neutron stars. 


1 INTRODUCTION 

Recently, a special attention has been paid to neutron star (NS)- 
NS or NS-black hole (BH) mergers (hereafter NSM), con¬ 
sidered today as the most promising sources of gravitational 
waves for ground-ba s ed detectors such a s advanced Virgo/ad 
Ligo jPh innevI 199 ll: Naravan et al. 1991 1 _ Wrgo^collaboratioM 


Abadie_et^ 'l2010l: LIGO Scientific c ollaboration I 201 ll : 


Aasi et al. I l2014al ll: d: iBelczvnski et al. 1 12014|) but also because 
several evidences (see lBergerll20l4 for a recent review) suggest 
in addition th at these events are the progenitors of short gamma-ray 
bursts (GRB) (|Pac zvnskill986l: lEichler et al. ll989l : lNaravan et aid 

ll992l : lMochkovitch et al. Il993h . Recent simulations show that, af¬ 

ter the merger, the relic BH-torus system can produce an ultra - 
relativistic ejection along its rotation axis dRezzolla et al. iFol ih . 
potentially leading to a short GRB. The similarities between the 
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properties of the pr ompt emission in short and long GRBs (se e.g. 
iGuiriec et al. suggest that th e same dissipation process is at 

work in these two classes of bursts dBosniak & Daign^l20l4) and 
that the main differences are the lifetime of the central engine, lead¬ 
ing in the case of binary compact objects to shorter events, and the 
density of the circumburst medium, leading to weaker afterglows 
dNakarlbOOTh . In addition, hydrodynamic simulations have con¬ 
firmed that a non-negligible amount of matter, typically about 10~^ 


to 10 ^ Mq, ca n be e j ected q uasi-isotropical l y i Janka_eLalJ[l999|; 


Rosswog^et^ 1^94 J2004; Oechslin et alj 2007 : Goriely et alj 
201ll: iBauswein et al.l 120131: iGorielv et al.ll2Q13l : I Just et alji2015F 


Wanaio et al.l2014l) and that the matter ejected in such NSM events 
is essentially made of r-process nuclei, shedding a new light on a 
old mystery. 

The r-process, or the rapid neutron-capture process, of stel¬ 
lar nucleosynthesis is invoked to explain the production of the sta¬ 
ble (and some long-lived radioactive) neutron-rich nuclides heav¬ 
ier than iron that are observed in stars of various metallicities, as 
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well as in the solar system (for a review, see lAmould et alj|2007h . 
The r-process remains the most complex nucleosynthetic process 
to model from the astrophysics as well as nuclear-physics points of 
view. The site(s) of the r-process is (are) not identified yet, all the 
proposed scenarios facing serious problems. Complex—and often 
exotic—sites have been considered in the hope of identifying as- 
trophysical conditions in which the production of neutrons is large 
enough to give rise to a successful r-process. 

Progress in the modelling of type-II supernovae and long 
GRBs has raised a lot of excitement about the so-called neutrino- 
driven wind environment. However, until now a successful r- 
process cannot be obtained ab initio without tuning the relevant 
parameters (neutron excess, entropy, expansion timescale) in a way 
that is not supported by the m ost sophisticated existing models 
dWanaio et alj201lHjankj|2012h . 


Early in the development of the theory of nucleosynthesis, 
an alternative to the r-pro cess in high-temperature su pernova en¬ 
vironments was proposed dTsuruta & Cameronll965h . It relies on 
the fact that at high densities (typically p > 10*° g/cm^) mat¬ 
ter tends to be composed of nuclei lying on the neutron-rich side 
of the valley of nuclear stability as a result of free-electron cap¬ 
tures. The astrophysical plausibility of this scenario in account¬ 
ing for the production of the r-nuclides has long been ques¬ 
tioned. It remained largely unexplored until the study of the de¬ 
compression of cold neutronised ma tter resulting from tidal ef¬ 
fects of a BH on a NS compan ion dLattimer & SchrammI 1 19741 : 
iLattimer et alJl9'7^ : lMevej|l989h . Many investigations with grow¬ 
ing sophistication have now confirmed NSM ejecta as viable 


2 OIII: IBauswein et al.ll2013l: iGorielv et alj 20131: Just et al.ll2015: 

Freiburghaus et alj 

I 999 I: Gorielv et al. 2005: Metzger et alJl201C : 

Roberts et al. 2011 

: Korobkin et al. 20121: Wanaio et al. 2014t). In 


show that the combined contributions of both the dynamical 
(prompt) ejecta expelled during the binary NS-NS or NS-BH 
merger and the neutrino and viscously driven outflows generated 
during the post-merger remnant evolution of the relic BH-torus sys¬ 
tems lead to the production of r-process elements from A > 90 up to 
thorium and uranium with an abundance distribution that reproduce 
extremely well the solar distribution, as well as the elemental distri¬ 
bution spectroscopically determined in very-low-metallicity stars. 
The ejected mass of r-process material, combined with the pre- 
dicted astrophysical ev ent rate (around 10“^ yr“* in the Milky Way 
jPominik et al. 2012h can account for the majority of r-material 


in our Galaxy jGorielv et al.ll201 ll : iBauswein et al.ll^l3l : llust et all 

[201^ . 


In NSM, nearly all of the ejecta are converted to r-process 
nuclei, whose radioactive decay heating leads to potentially ob- 
servabl e electromagnetic radiation in the optic al and infrared 
bands jLi & Paczviiski|[l9^ : iMetzger et al.ll2010l) with 100-1000 
times fainter peak brightn esses than those of typical supemovae 
and durations of only days l lGorielv et al.l20li] : lRoberts et alj2oTII : 
IBauswein et al.|[^13tl . These “macronovae” or “kilonovae” are in¬ 
tensely searched for (with a recent, possible first success in as- 
sociation with the short GRB 130603B, see iBerger et al. lEoiSl : 
iTanvir et al.ll2013l) . Their unambiguous discovery would constitute 
the first detection of r-material in situ. 

In the present paper, our primary goal is to predict the cos¬ 
mic evolution of the NSM that can be deduced from the ob¬ 
served history of the star formation rate (S FR), including the re¬ 
cent constraints obtained at high redshift JBehroozi_et^ 201G; 
iBehroozi & SilkI l2015l : lOesch et"^ l2014l : iBouwens et al I 2014 ; 


iKistler et al.ll2013l : iTrenti et^l2013l) . Computing the NSM rate 
requires to determine poorly known parameters, namely the frac¬ 
tion of NS in a binary system with another NS or a BH, and the 
characteristic coalescence timescale of such a system. To constrain 
these parameters, we model the associated evolution of a typical 
r-process element. Europium, by assuming NSM as the main astro- 
physical site of the r-pr ocess and adopting the most recent nucle- 
osynthetic calcu lations jGorielv et aljboill : IBauswein et alj|2013l : 
llust et al1l2015l) . and we compare to the available observations in 
old stars of the Galactic halo and in external galaxies. Once the 
merger rate has been determined, it is possible to estimate the ex¬ 
pected event rate for the new generation of gravitational wave de¬ 
tectors (advanced Virgo and LIGO), as well as the expected de¬ 
tection rate of kilonovae in present and future optical/near-infrared 
large surveys. The same approach can be followed to determine, 
without any new parameter, the CCSN rate, hence the predicted 
evolution of Europium in the alternative scenario where r-process 
elements are produced mainly in such stellar explosions. 

The paper is organized as follows: in Sect. |2l we select three 
possible scenarios for the SFR, which are representative of the cur¬ 
rent uncertainties at high redshift, and we describe how the CCSN 
and NSM rates can be deduced from the SFR, both scenarios still 
depending on given unknown parameters. In Sect[3 we describe 
the chemical evolution model that we use to deduce from the cos¬ 
mic SFR the expected evolution of the abundances of various chem¬ 
ical elements in star-forming structures. We show that the three sce¬ 
narios considered in the paper for the cosmic SFR are compatible 
with a large set of observations, hence validating the present model. 
In § 14.21 we show how the model can be extended to heavy elements 
produced by the r-process. In Sect.|4l we present the predicted cos¬ 
mic evolution of Eu, assuming either the CCSN or the NSM as the 
main astrophysical site of the r-process. We make a detailed com¬ 
parison with observations at low metallicities which correspond to 
the early evolution in the Universe and appear to be the most dis¬ 
criminative between both sites. In the NSM case, this allows us in 
addition to constrain the properties of the binary system, in partic¬ 
ular coalescence timescale. The resulting lower and upper limits on 
the merger rate are used in Sect[5]to make predictions for the ex¬ 
pected event rate in the horizon of advanced Virgo/LIGO, as well 
as for the kilonova rate in future surveys. Finally, we draw our con¬ 
clusions in Sect.[^ 


2 PREDICTING THE NSM RATE 
2.1 The cosmic star formation rate (SFR) 


The cosmic SFR is a key ingredient to model the rate of stellar ex¬ 
plosive events, such as CCSN and NSM. The redshift evolution of 
the SFR density is now rather well constrained by many observa¬ 
tions. Recent data from the Hubble Ultra Deep Field have signif¬ 
icantly extended the ra nge of redshift for its determination, from 
redshift z = 4 up to 10 1 Boi^ens^t a ^20071 : iBouwens et al 11200^ 


I 2 OI lLl20l4l : IOesch et al.ll2012l . l2013ll2014l) . Even more recent ob¬ 

servations of high z galaxies and GRBs tend to favor a large a mount 
of still unseen SFR at z > 9 jKistler et al.ll2013l : IWandl2()13l) . This 
last study suggests that SFR density may only decline out at z = 11 
and that GRBs may be useful in explor i ng the unseen faint dwarf 
galaxies at high redshift. [Behroozi et alj l l2013tl re-analyzed the av¬ 
erage star formation histories from z = 0 to 8 and obtained con¬ 
sistent results with observed galaxy stellar mass functions, specific 
star formation rates and SFR. Moreover, high-redshift galaxy evo- 
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Time (Gyr) 

13.67 5 1 0.5 



For each mode, the IMF slope is set to the Salpeter value, i.e. x = 
0.1024 1.35. 

Star formation is assumed to start at the initial time of 1q = 
100 Myr corresponding to a redshift of z = 30. We adopt the 
following cosmological parameters, 0,„ = 0.27, Ha = 0.73 and 
Hq - 1\ km/s/Mpc (h = 0.71) and a primordial power spectrum 
with a power-law index « = 1. The age t of the Universe is then 
related to the redshift hy 

dt 9.78 Gyr 

— = - - - ( 2 ) 

dz (l+z)[nA-m,(l-tz)3]“'5 

As seen in Fig. [T] the SFR2 mode is an intermediate case, 
whereas SFRl and SFR3 are respectively close to the lower and 
higher values of the current measurements at high redshift. We 
checked that each of these three scenarios is compatible with a large 
set of independent observational constraints discussed in Sect. l3.2l 
especially the Thomson optical depth to the cosmological micro- 
wave background (CMB) and the cosmic evolution of [Fe/H]. This 
led to the inclusion of the additional Pop III component at high red¬ 
shift in the SFRl case, the standard mode alone being too weak in 
this case to fulfill the reionization constraint (see Sect. l3.2.7T l. 


Figure 1. Cosmic SFR as a function of redshift. Three SFR modes consid¬ 
ered in the paper, SFRl (solid line), SFR2 (dotted line) and SFR3 (dashed 
line). Observat i ons a r e taken from iBehro ozi et al.l l2013h (red points), 
ISouwens et^ 1 j2014h:IOesch etalJ <2014 (and references therein) (blue 
Dointsl and iKistler et alJ<20l4 (black points). 


lution has also been predi cted by Behroozi_^Sil^ 2015h who, in¬ 
cluding recent constraints dSehroo zi et alJ2013l : IOesch et al.l201.3l : 
iKistler et alj2013l : lTrenti et alj2013l) . probed~the unobserved z > 8 

galaxy populations. The resulting observationally constrained SFR 
is illustrated in Fig. [T] as a function of the redshift. The corre¬ 
sponding evolution of the cosmic SFR density as a function of 
redshift can be parametrize d by the following form proposed by 
ISoringel & HemauistI ll200.3h . 


ip(z) = V- 
a 


flexp[fc(z-z„,)] 

-b + hexp[a(z - z™)] ’ 


( 1 ) 


where v (in Mo/yr/Mpc^) and Zm correspond to the astration rate and 
the redshift at the SFR maximum, respectively, while b and b - a 
fix the SFR slope at low and high redshifts, respectively. 

Due to the fundamental importance of the SFR in predicting 
the stellar explosive event rates, and more generally the cosmic 
chemical evolution, we consider in the present paper three SFR 
modes (see Fig. [T}, which coincide at low redshift, where they re¬ 
produce the available data, and differ at high redshift, taking into 
account the current uncertainties: 


• SFRl includes a standard mode of population (Pop) II/I stars 
formation between 0.1 and 100 Mo, and additionally a Pop III stel¬ 
lar mode at high redshift between 36 and 100 Mq. The correspond¬ 
ing parameters for each mode in Eq. 0 are, respectively, v = 0.18 
and 0.0025, z,,, = 2.0 and 12.0, a = 2.37 and 4.0, and b = 1.80 and 
3.36. 

• SFR2 includes a unique mode of star formation between 0.1 
Mo and 100 Ma, with the following parameters: v = 0.15, Zm = 1-7, 
a = 2.8, b = 2.45; 

• SFR3 also includes a unique mode of star formation between 
0.1 Mo and 100 Mq, with the following parameters: v = 0.36, z,,, = 
2.6, a = 1.92, h= 1.5. 


2.2 The NSM rate 

2.2.1 The CCSN rate and the NS birth rate 


Once the evolution of cosmic star formation rate is known, the 
CCSN rate can be directly deduced if the ma ss-range of the pro- 

f enito rs is known. We use the predictions from lWooslev & Weaver! 
1995h and assume that the remnant of stars with mass 8 < M/Mo < 
30 is a NS, and that more massive stars produce a BH. This leads to 
the following NS birth rate at time t (related to redshift z by Eq.0: 

RNs(r)= {Amtb{mmQ^ns{m,Z{t.)) , (3) 


where is the initial mass function (IMF) integrated over the 
whole stellar mass range, T(r) the SFR, t, the formation time of the 
progenitor, such that t,+T (m, Z(t,)) = t, Z{t,) the metallicity at the 
epoch of formation of the progenitor, t (m, Z(t,)) the lifetime of a 
star formed with mass m and metallicity Z(f,), and Hns (ni, Z(t,)) is 
1 if the stell ar remnant is a NS, 0 othe rwise. The lifetimes of stars 
come from jMaeder & Mevnelll989l) and the stellar metallicity - 
dependent remnants are taken from dWooslev & Weav^ Il995l) . 
The BH birth rate is given by a similar equation, using now 
3 bh (ttt, Z(t,)) = 1 if the remnant is a BH, 0 otherwise, and the 
CCSN rate is the sum of these two rates. The corresponding evo¬ 
lution of these three rates deduced on the basis of SFRl and SFR2 
is plotted in Fig. |2] On the observational point of view, the local 
CCSN rate is well determi ned, thanks to large surveys dedicated 
to the search of supernovae dMattila et al.ll2012h . As can be seen in 
Fig .|2l the agreement with the prediction deduced from the cosmic 
SFR is good. 


2.2.2 The NS merger rate 

In contrast to the to the CCSN rate, the astrophysical esti¬ 
mate of compact-binary coalescence rates depends on a num¬ 
ber of assumptio ns and unknown mo del parameters and are still 
rather uncertain dAbadie et al.l l2010l) . Among the various esti¬ 
mates, the rate predictions for NS-NS binary systems is the most 
reliable one since they are based on extrapolations from ob¬ 
served binary pulsars in our Galaxy. These estimates lead to 
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IPominik et alj (1201 3h calculated the cosmological merger rates of 
NS-NS, NS-BH, BH-BH systems as a function of redshift assum¬ 
ing different SFR histories. While in most cases NS-NS systems 
dominate the merger rates in the local Universe, BH-BH mergers 
are found to dominate at high redshift. 

In the present study, the merger rate is determined assuming 
that the formation rate of NS-NS/BH binary systems is a fraction a 
of the NS birth rate Rns and that these binary systems merge after 
a delay At^sM. the value of which is discussed below. This leads to 
the following merger rate at time f. 

^nsm( 0 = a^NS (t “ AfNSM) • (4) 

In our model, both a and A/nsm and taken as free parameters. The 
predicted evolution of the merger rate is plotted in Fig.j^for SFRl 
and SFR2, assuming a = 0.002 and AJmsm = 0.2 Gyr. The inde- 
pendent estimate o f the local rate taken from the compilation by 
lAbadie et alj (1201 Oh is indicated for comparison. 


2.2.3 The coalescence timescale 



Redshift z 


a, 



S 


K 

pH. 

tZ3 


Figure 2. NS and BH birth rates as well as the CCSN and merger rates as 
a function of the redshift for SFRl (upper panel) and SFR2 (lower panel) . 
Green and red points represen t the local CCSN r ate jMattila et alj|20l3) 
and the total local merger rate iAbadie et al.ll201oli . respectively. The SFR 
evolution is also plotted for comparison. 


a coalescence rate of about 100 Myr^' per Galaxy, although 
this rate could plausibly range from 0.1 Myr^' to 1000 Myr^* 
dKim et alJl2003l : iKalogera et alJl2003 : iMennekens & VanbeverenI 


12014 ). Due to the lack of observational data, predictions of the 
frequency of NS-BH mergers remai ns even more uncertain and 
rely exclusively on theoretical studies (|Tutukoy_&_VungeNoiiLL993; 
Voss & T aurij l20 03l: lO’Shaughnessv et al.l 1200^: iDom inik et al l 


201 2 I: Dominik e t al J~ 2013 : Mennekens & VanbeverenI 20141 : 

Postnov & YungelsonI 1201 jk These population synthesis models 


estimate galactic merger rates between 2 x 10“^ and 10“^ per year. 
The range reflects the challenge in comprehensively modelling 
the formation and evo lution of stellar binaries and their remnants. 
iBauswein et al.l ( 12014l) determined an independent upper limit on 
the merger rate of NS-BH binaries by comparing the predicted r- 
process nucleosynthesis yields of such systems with the observed 
galactic amount of r-process material. A strict upper limit of the 
average NS-BH merger rate of about 6x10“^ per year was found. 


The coalescence timescale AJnsm is defined as the delay between 
the formation of the binary system of two NS or a NS and a BH 
(note that it does not include the lifetime of the two star progeni¬ 
tors of the compact objects). It is a key ingredient as it leads to a 
different cosmic evolution of NSM compared to CCSN, for which 
there is no additional delay after the formation of a NS/BH. 

The angular momentum loss due to the emission of gravita¬ 
tional waves governs the coalescence timescale. Its value strongly 
depends on the properties of the binary system, and in particular on 
the initial separ ation a (or orbital period T ) since Af ^sM (resp. 
AtNSM itPsters & Mathews lll963L see also IKalogera et akl 

l200ll : lHughes II 2 OO 9 I) . This dependence leads to an expected broad 
range of values for At^sM. with however large uncertainties related 
to the distribution of pos sible initial orb ital parameters in NS/NS 
and NS/BH binaries (see iHughes II2009I . for a recent review). For 
a binary system of two 1.4 Mq NS, the coalescence timescale is 
AfNSM = 64 Myr for a = 10“^ AU (resp. T = 5.2 h). Only 
seven NS/NS binary systems are known with measured masses 
dUorimer I2005ll200^ . The current orbital parameters and the mea¬ 
sured masses allow us to estimate the remaining time before the 
merger. As at least one of the two NS is observed as a pulsar, 
the measurement of the pulsar period and its time derivative can 
provide an estimate of the current age of the binary system. The 
sum of these two t imescales gives an estimate of AfNSM- Using 
the values listed in iLorimer I J2008l) . coalescence timescales be¬ 
tween 100 and 400 Myr are found in four cases, and values larger 
than 1 Gyr in the three other cases. The double pulsar PSR 1037- 
3039 corresponds to the lowest value, Af^sM - 180 Myr. The fa¬ 
mous binary pulsar PSR B1913-I-16, which has been studied in de¬ 
tails for years as a test of genera l relativity dHulse & Taylor 1^9751 : 
IWeiberg, Nice and Taylor IItOIOI) is characterised by a coalescence 
timescale AJnsm - 420 Myr. 

The distribution of the coalescence timescale can be esti¬ 
mated using a population synthesis code for binary systems, which 
can be calibrated with observed systems. The result is plotted for 
two different metallicities in Fig . [3l It has been obtained usin g 
the StarTrack code described in l^lczvnski et all d2002Ll^08al) . 
where standard stellar evolution models are implemented with the 
most updated physics regarding the relevant processes for the evo¬ 
lution of binary systems. Fig. [^has been produced using the re¬ 
sults of a recently updated version of the code (Mink & Belczyn¬ 
ski in preparation), including new prescriptions for stellar wind 
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Figure 3. Distribution of the coalescence timescale A^nsm for NS/NS 
binary systems. The coalescence timescale A/nsm is defined as the delay 
between the formation of the two NS and the final merger. The plot shows 
the fraction of mergers with a coalescence timescale above a given value 
in Myr, for two different metallicities. This distribution has been obtained 
from an updated versio n of the population synthesis model described in 
iBelczvnski et alJ j2002h . It covers a broad range of values with a distribution 
falling approximately like p(AtNSM) ^^nsM' 


mass loss rates l[^nketa^ 20n; Vink & Grafenerll2Q12h . for su- 
pernova models jFrver et al, 12012 ), for the common envelop phase 
toominik et al. 1120 12 ), adding elect ron capture supernovae for NS 
formation teelczvnski et al]|2008ah . and taking into account re- 
cent observationa l constraints on massive stars in binary systems 
jSana etal. ,l2012h . Additionally, the evolutionary model included in 
the code does not allow for the progenitors to undergo a common- 
envelope phase while one of the massive binary stars is evolving 
through the Hertzsprung gap fsee lBelczvnski et al. II2007I . for more 
details). This corresponds to a conservative assumption and ad¬ 
ditional NS-NS mergers m ay form if this assumpt ion is relaxed. 
Compared to the results of IBelczvnski et aP ( I 2 OO 2 I) . this updated 
calculation predicts a fraction of undetected galactic merging NS- 
NS systems significantly smaller, now at the level of ~ 10 %. As 
illustrated in Fig. [3] 

the distribution of coalescence timescales predicted by such a 
detailed calculation shows a broad range of values: at solar metal- 
licity Z = Zq (resp. at metallicity Z = 0. IZq), it is found that 
29% (resp. 28%) of systems have AfNSM < 100 Myr, 13% ( resp. 
16%) have 100 < ArNSM ^ 300 Myr, 15% (resp. 14%) have 
0.3 < AfNSM ^ 1 Gyr, 17% (resp. 16%) have 1 < AfNSM ^ 3 Gyr, 
and 26% (resp. 26%) have AfNSM > 3 Gyr. The impact of the metal¬ 
licity is small. It appears clearly that a significant fraction of sys¬ 
tems has a coalescence timescale of a few 100 Myr or less. These 
systems may dominate the early evolution of the NS merger, and 
therefore the evolution of Eu at redshift discussed in Sect.O 


3 COSMIC CHEMICAL EVOLUTION 


3.1 Model description 


In addition of the rates of explosive events (CCSN, NSM) discussed 
in the previous sections, the predicted chemical evolution can also 
be deduced from the cosmic evolution of the SFR. We use the 


model develop e d by Daigne et al. I j2006l) : iRollinde et all ( l2009l) : 
IVangioni et alj (12014 1. which is based on a hierarchical model 


for structure formation i Press & SchechteJ 

1974:ISheth & TormenI 

Il999l:ljenkins et al. 11200ihlWvithe & Loeb 

|2003f). We assume that 


the minimum mass of dark matter haloes for star-forming structures 
is 10 ’ Mq. 


Finally, for a given evolution of the cosmic SFR density, as 
discussed in § 12.11 the model follows the evolution of baryons in 
the Universe in terms of two main reservoirs. The first is associ¬ 
ated with collapsed structures and is divided in two sub-reservoirs: 
the ISM gas (mass MisM(f)) and the stars and their remnants (mass 
The second reservoir, the intergalactic medium (IGM, mass 
AIiGM(f)), corresponds to the medium in between the collapsed 
structures. The evolution of the baryonic mass of these reservoirs 
is governed by a set of differential equations: 


dMiou 

dt 


dl^atiuct 

dt 


-flbfr) -I- o(f). 


(5) 


dM, 

dt 


- 4^(0 - e(f) and 


dt 


dM,trua _ dM, 
dt dt 


(6) 


In addition, we have Mism(0 + M,{t) = MstmctfO. corresponding to 
the total baryonic mass of the structures, and MiGuit) + Afstruct(f) = 
constant (total baryonic mass of the Universe). 

As can be seen, these equations are controlled by four rates 
which represent four fundamental processes: the formation of stars 
through the transfer of baryons from the ISM, Tfr), which is di¬ 
rectly fixed by the assumption SFRl, 2 or 3 discussed in § 12.11 
the formation of structures through the accretion of baryons from 
the IGM, flb(0; the ejection of enriched gas by stars, e{t) and the 
outflow of baryons from the structures into the IGM, o(t). 

In the model, we track the chemical composition of the ISM 
and the IGM as a function of time or redshift. The differential equa¬ 
tions governing the evolution of the mass fraction (A'°t^) of 
element i i n the ISM (IGM) are given by equations (6) and (7) in 
section 2 in lOaigne et al. I l l2004ll . 

The baryon accretion rate, flb(f) is computed in the framework 
of the hierarchical scenario of structure formation and the baryon 
outflow rate, o(f), from the structures includes a redshift-dependent 
efficiency, which accounts for the increasing escape velocity of the 
structures as the galaxy assembly is in progress. 

We use the framework of the hierarchical scenario where small 
structures are formed first. At redshift z, the comoving density of 
dark matter halos in the mass range [M, M + dM} is fps{M, z)dM, 
with 



dM M fps(M,z) - Pdm . 


(7) 


where puM is the comoving dark matter density. The distribution 
fu nction of halos fps(M , z) is computed using the method described 
in ijenkins et al. I l l200lh . We adopt a rms amplitude trg = 0.9 for 
mass density fluctuations in a sphere of radius 8 Mpc. We as¬ 
sume that the baryon distribution traces the dark matter distribution 
without any bias so that the density of baryons is just proportional 
to the density of dark matter by a factor n*/ (Q.„, - Q/,). 

The fraction of baryons at redshift z which are in such struc¬ 
tures is given by 
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r dMMfys(M,z) 

_p 

/b,struct ~ roa ■ 

dMMfps{M,z) 

Therefore, the mass flux can be estimated by 



/dtr' 

struct 


IdzJ 

dz 


( 8 ) 


(9) 


In addition to the NS and BH birth rates, the CCSN rate, 
and the NS-NS and NS-BH merger rates discussed in § 12.21 land 
of special interest in the present study), several other quantities 
are followed as a function of the redshift z, namely the WD birth 
rate, the type la supernova rate, the abundances of various chem¬ 
ical elements, including Fe and Eu, in the ISM and the IGM, 
the ionizing flux from stars, the ionization state of the IGM and 
the Thomson optical depth of the CMB. A detailed description o f 
this cosmic evolution model can be found inlPaigne et al. UTOObh : 
iRollinde et alJ (l2009h : IVangioni et aTI i2014h . 


3.2 Observational constraints on the chemical evolution 
model 

We tested that each cosmic evolution scenario discussed in this pa¬ 
per fulfills two important observational constraints, namely the ob¬ 
served cosmic evolution of [Fe/H] (corresponds to the logarithm 
of the abundance of iron normalized to the solar one) and the 
measured Thomson optica l depth of the CMB. Following the de¬ 
tailed analysis presented in IVangioni et al.l ll2014h for different pos¬ 
sible SFR using the same cosmic evolution model, this leads to the 
choice of the three SFRs presented in § 12.11 The SFRl and SFR2 
cases are directly taken from their study. 

SFRl corresponds to a standard SFR with an additional Pop 
III star component at high redshift, mainly motivated by the fact 
that the ionizing flux at high redshift produced by the standard 
mode only is not large enough to reproduce the Thomson optical 
depth of the CMB (see § [Tm (. 

SFR2 is a SFR enhanced at high z. The third mode considered 
here, SFR3, is chosen as the upper limit of the SFR at high redshift 
obtained from the constraint of the Thomson optical depth of the 
CMB (see Fig. [^. As shown in Fig.[Tl SFRl fits the observat ional 
constraints fromiBehroozi et alJ ll2013h : lBouwens et al I ll2014h and 
lOesch et al.U2014h. whereas SFR2 and SFR3 take into account the 
studies of jTrenti et alj2013l) and iKistler et alj ( l2013h . respectively. 


3.2.1 Reionization 


Concerning the ionization history, the evolution of the volume fill¬ 
ing fraction Qmniz) of ionized regions at redshift z is given by 


d2ion(z) 1 d/iic 


dz 


nb dz 


—-apn^C{z)QL{z){l+z? 


( 10 ) 


where nj, is the comoving density in baryons, nion(z) the comov¬ 
ing density of ionizing photons, ap the recombination coeffi- 
cient, and C(z) the clu mping factor. This last factor is taken from 
ICreif & BrommI ( l2006h and varies from a value of 2 at z < 20 to a 
constant value of 10 for z < 6. The escape fraction, /esc, is set to 0.2 
for each of our SFR modes. The number of ionizing phot ons emit¬ 
ted by massive stars is calculated using the tables given in dSchaereJ 
1 20021) . Finally, the Thomson optical depth between redshifts z' = 0 
and z is computed following dCreif & Brommll200(il) . i.e. 


T - CCr-p H'a 


r dz' 

Jo 


QUz')(l+zf 


(z') 


( 11 ) 



0 5 10 15 20 

z 



Figure 4. Reionisation. Evolution as a function of redshift of the ionizing 
flux (top panel), volume filling fraction (middle panel) and Thomson 
optical depth (lower panel) for the three SFR modes (conventions are the 
same as in Fig. [T). W MAP9 results for the optical depth at high redshift 
iFlinshaw et alj201^ are indicated by a red strip in the lower panel. 


where crp is the Thomson scattering cross section. Fig. |4] shows 
the evolution of the volume filling fraction Qion and the ionizing 
flux for the three SFR histories, and the resulting Thomson optical 
depth, the value of which is found to agree well with the measure- 
ment from CMB observations by WMAP9 dHinshaw et alJlTOllh 
at high redshift. As mentioned above, the SFRl scenario without 
the Pop III mode at h igh redshift would not fulfill this constraint 
dVangioni et al.ll2oT^ and SFR3 is adjusted to reproduce the up- 
per limit on the optical dep t h. No te that the last results from the 
IPlanck Collaboration I d20l4l20T5h give a revised Thomson optical 
depth T = 0.079 ± 0.017. This value is lower than the WMAP9 one 
( T = 0.089 ± 0.014). Consequently, in this context, the contribution 
of Pop III stars in the SFRl model should be lowered. 
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3.2.2 Cosmic chemical evolution: iron 

Our model incl u des t he stellar lifetimes estimated by 
iMaeder & MevneJ ( Il989ll for interme diate mass stars with 
0.9 < M/Mo < 8 and bv ischaererl l l2002h for more massive stars. 
Lifetimes of low mass stars (M/Mq < 0.9), such as the old halo 
stars where abundances of Eu are measured at low metallicity, 
are long enough for the star to be still observed today, whatever 
the redshift of formation. Stars formed at redshift z inherit their 
initial chemical composition from the abundances in the ISM at 
the same redshift. Thus, the observed abundances in stars reflect, 
in a somewhat complex way, the yields of all massive stars that 
have exploded earlier. The stellar yields adopted in our model 
depend on the stellar ma ss and metallicity. By de fault, we use 
the tables of yields from dWooslev & Weav^ll995h for massive 
stars (8 < M/Mq < 40). An interpolation is made between the 
different metallicities (Z = 0, O.OOOIZq, O.OOIZq, O.IZq and Zq) 
and masses, and tabulated values are extrapolated for stars beyond 
40 Mq, which correspond to a small fraction of the population of 
massive stars (i.e. (40“-' - 100“') / (8“-' - lOO^-') 8% for the 

Salpeter IMF of index x = 1.35). 

Fig. [5] shows the [Fe/H] evolution as a function of time and 
redshift for the three SFR modes. Data correspond to chemical 
abun dance determinations in about 250 damped Lya (DLA) sys¬ 
tems jRafelski et al.ll20O . and references therein) as a function of 
the redshift. In their paper, they present 47 DLA systems at z > 4 
observed with the Echellette Spectrograph and Imager and the High 
Resolutio n Echelle Spectromete r on the Keck telescopes (for more 
details see lVangioni et alj2014h . We use Eq. (|3 to associate a time 
with the observed redshift for each observed DLA system. As il¬ 
lustrated in Fig. [3 the cosmic evolution of the iron abundance is 
reproduced by each of the three scenarios, passing through the up¬ 
per limit of the points. However, note that these iron abundances 
are measured in the gas phase which should be considered as lower 
limits due to depletion into grains in the ISM. As seen in Fig. |5] 
(middle panel), the early iron evolution depends strongly on the 
shape of the SFR at high z. In particular, it can be seen that a cos¬ 
mic [Fe/H] = -4 enrichment is achieved at the very beginning, 
At = t-to=^ 0.015 Gyr for SFR2, 0.05 Gyr for SFR3, and 0.4 Gyr 
for SFRl. This time delay in the early Fe enrichment can have a 
drastic impact on the interpretation of the early cosmic chemical 
evolution with respect to observations, as discussed in Sect.|4] 

3.2.3 Cosmic chemical evolution: other elements 

As detailed in IVangioni et alj (l2014h . the present cosmic evolution¬ 
ary model has been used to predict the “standard” chemical evolu¬ 
tion of light elements, such as C, N, O or Mg and to test its predic¬ 
tive p ower with respect to available observations jSuda et al.l2008L 
l201lh . The sensitivity to different SFR has also been explored. 

• Carbon, oxygen. It was found that the [C/H] evolution as a 
function of the iron abundance is essentially independent of the 
adopted SFR model for [Fe/H] > -2 or z < 4. At higher z (lower 
metallicity), the SFR including a Pop III mode (like SFRl) was 
shown to be in better agreement with the observation of ultra- 
metal-poor carbon-enhanced stars, which presents a significant dis¬ 
persion. Similar conclusions were drawn for the oxygen evolution, 
where the Pop III contribution to the SFRl model allows for a high 
production of O at low metallicities with a local peak in [O/H] 
at [Fe/H] sj -4.8. Interestingly, this model can also explain the 
oxygen abundances in ultra-l ow-metallicity stars (see in particular 
Fig. 9 in IVangioni et alj2014l) . 
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Figure 5. Iron evolution. [Fe/H] evolution as a function of time t (upper 
panels) or redshift z (lower panel) for the three SFR modes (lines are the 
same as defined in Fig. \Q. The middle panel focusses on the early time 
evolution, with the insert emphasizing the thi'ee different time delays around 
the [Fe/H]=-4 iron en richment . Data corr espond to DLAs observations as 
a function of redshift ^Rafelski et al.l2Ql3) . Eq. (2} is used to convert z into 
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• Magnesium. As in the oxygen case, the inclusion of a high- 
mass model in SFRl at large redshift directly impacts the evolution 
of Mg in star forming structures. The predicted Mg pattern in the 
SFRl case is, hovyever, no t observed at low metallicity (see Fig. 
10 in I Vangioni et al.l|2014h . However, new studies in high mass 
stars show that Mg evolution at low metallicity remains sensitive 
to the predicted nucleosynthesis in massive zer o-metallicity stars. 
Such c alculations have been recently revisited bv lHeeer & Woosl^ 
( l2010h who showed that in their models no Mg is produced in 
zero-metallicity stars more massive than 30 Mq. In this case, the 
early production of Mg in the three SFR cases discussed here is 
found to reproduce w ell the bulk of data, as shown in Figure 10b of 
IVangioni et alj ( l2014h . 

• Local stellar metallicity distribution fun ction (MDF). 

We show in Fig. [6] the MDF derived by lAn et alj ( 1201 3h us¬ 
ing the SDSS photometric Catalogue together with our computed 
MDF based on the SFRl m ode (the other modes giving very 
similar results). As found by lAn etalJ i l2013ll . the in-situ photo¬ 
metric metallicity distribution has a shape tha t matches that of 
the ki nematically-selected local halo stars from lNorris and Ryan I 
( Il99ih . For this reason, we adopt their SDSS calibration catalog 
MDF as a reference since it presents the advantage of expressing 
the MDF directly as a function of [Fe/H]. The global shape of our 
computed MDF is shown in Fig.[^to be in good agreement with the 
observation. We point out here the impact of the choice of the iron 
yields and the related uncertainties. We present four curves reflect¬ 
ing uncertainties mainly related to the CCSN mass cut. The blue 
curve corresponds to the standard yields from lWooslev & Weaved 
( Il995h adopted in the rest of the paper, which are varying with 
metallicity, the green curve to the same yields divided by a factor 
2 , the dot-dash black curve to yields calculated from solar metal- 
licity stars only , and the red curve to Fe yields determined by 
iKobavashi et ^( l2006h . 

• Stellar mass. As shown bv IVangioni et alj l l2014h . the cosmic 
evolution of the stellar mass as a function of the redshift is also 
found to be in good agreement with cosmological observational 
photometric data (see their Fig. 12b), in particular for the SFRl 
mode. 

_ In conclusion, based on the previous study of IVangioni et all 

( l2014h within the same cosmological context and on the basis of 
the same SFR cases, it can be concluded that the standard chem¬ 
ical evolution is rather well reproduced by the cosmic evolution 
of the star formation considered in the present study, and conse¬ 
quently that the present model is well suited to describe the av¬ 
eraged chemical history. The evolution of Milky Way-like galax¬ 
ies appears to be not far from the a veraged history. As discussed 
recently in iFrebel and Norrisl l l2015l) . this opens the possibility of 
"near-held cosmology” based on stellar archeology, especially us¬ 
ing the population of the oldest stars in the halo of our Galaxy and 
in dwarf galaxies. We will use this approach in the next Section to 
study the cosmic evolution of r-process elements. 


4 COSMIC EVOLUTION OF R-PROCESS ELEMENTS 

As far as the galactic chemical evolution is concerned, CCSN re¬ 
main promising sites, especially in view of thei r potential to signif - 
icantly contribute to the galactic enrichment dArgast et al.ll2004l) . 
However, they remain handicapped by large uncertainties associ¬ 
ated mainly with the still incompletely understood mechanism that 
is responsible for the supernova explosion and the persistent dif- 
hculties to obtain suitable r-process conditions in self-consistent 



[Fe/H] 


Figure 6. Local metallicity distribution function, i.e. normalized number 
of stars as a function of the iron abundance [Fe/H]. The MDF is calcu- 
late d with SFRl m ode and the red circles to the observed MDF deduced 
by dAn et alj|2013h using SDSS photometry. We present four curves ob¬ 
tained with different Fe yields reflecting uncertaintie s mainly related to 
the C CSN mass cut. Blue curve : standard yields from dWooslev & Weav^ 
I 1995 I 1 with yields varying with metallicity, green curve : same after divid¬ 
ing these yields by a factor 2, dot-dash black curve : same using only yields 
calculated from solar meta llicity stars only, red cu rve : same when using the 
iron yields determined by dKobavashi et alj2006ll . 


dynamical explosion and neutron-star coo ling models ll.lankal 20121 : 
iHtidenohl et al.llToi(]l : lFischer et al.ll201C)h . In addition, nucleosyn¬ 
thesis predictions of the detailed composition of the ejected matter 
remain difficult due to the remarkable sensitivity of r-process cal¬ 
culation to the still unknown initial properties of the ejecta. 

In contrast, NS-NS or NS-BH mergers, which today are 
clearly favored from a nucleosynthesis point of view, have been 
claimed to be ruled out as the dominant r-process source by early 
stu dies on the basis o f inhomogeneous chemical evolution mod¬ 
els dArgast et alj|2004) . due to their low rates of occurrence which 
seem to be inconsistent with observations of low-metallicity r- 
process-rich stars. In addition, the significant amount of r-process 
material ejected by a single NSM leads to a large scatter in r- 
process enrichment at later times that does not seem to be con¬ 
firmed by observations. _ _ _ 

R ecent studies 1 Matteucci_et ^ |20 i 4 jKomi^^e^ffi 


20^; lCescutt^&^hiagginiJ_j201^]_jMennetens_^^nbeverer 


2014 Ir suiimoto & Shire vamal I2014alf9: IShen et alj 201 4 
de Voort et al. 20151 : iRamirez-Ruiz et al.l 1201 Sh h^e recon¬ 


sidered the chemical evolution of r-process elements in different 
evolutionary cont exts, and reached rath er different conclusions. 
More specifically. iMatteucci et al] d20l4) explored the Eu produc¬ 
tion in the Milky Way using a local chemical evolution model. 
Their observational constraints at low metallicity c ome essentially 
from t he obseryations of Eu in metal-poor stars by [Francois et al.l 
d2007l) . without including the related error bars and upper limits. 
The relevance of the NSM scenario on the production of Eu has 
been studied by testing the effect of (i) the coalescence timescale 
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of the binary system (ii) the Eu yield expected to be ejected from 
NSM and (Hi) the range of the initial mass of the NS progenitors. 
Similarly, the CCSN scenario has been explored by considering 
different possible Eu yields assuming a strong r-process taking 
place in CCSN. In this framework, NSM is found to be potentially 
a major r-process source if the coalescence timescale is of the 
order of 1 Myr and the ejected Eu yield of the order of 3 x 10“^ 
Mo for a mass range of progenitors of NS ranging between 9 and 
50 Mq. The scenario where both CCSN and NSM contribute to 
the Eu synthesis is also compatible with observations provided 
NSM produce 2 x 10“’ Mq of Eu per system and each CCSN with 
progenitors in the range of 20-50 Mq produce around 10“* - 10“^ 


Mo of Eu. 

In parallel. [Komiva et alj ( l20l4) investigated the chemical en¬ 
richment of r-process elements using a hierarchical galaxy forma¬ 
tion model. The CCSN scenario is found to reproduce the scatter 
of observed r abundances in low metallicity stars if about 10% of 
CCSN in the low-mass end (i.e. for progenitor mass of the order 
of 10 Mo) is the dominant r-process source and the star formation 
efficiency amounts to about 0.1 per Gyr. Eor NSM to be the main r- 
process site, a coalescence timescale of about 10 Myr with an event 
rate about 100 times larger than currently observed in the Galaxy 
need to be considered. 

We note t h at the coal escence timescale in the studies by 
iMatteucci et al.l ( l2014h and iKomiva et alj j2014l) is constrained 
to be surprisingly short compared to known systems in our 
Galaxy and theoretical predictions by stellar population models 
(see Sect. 12.2.31. _ 

ICescutti & Chiappini I l l2014l) computed inhomogeneous 
chemical evolution models for the galactic halo, taking into ac¬ 
count the contribution of electron-capture and magnetorotationally 
driven supernovae (including fast-rotating progenitors), but not 
NS M, to explain the Eu scatter i n meta l-poor stars. 

iMennekens & VanbeverenI ( l2014h also studied the temporal 
evolution of the galactic population of double NS binaries, mixed 
systems with a NS and BH component, and double BH binaries. 
They conclude that, except for the first 100 Myr of the evolutionary 
phase of the Galaxy, double compact star mergers may be the ma¬ 
jor production sites of r-process elements, and it is probable that the 
mixed NS-B H systems domina te o ver double NS binary merg ers. 

Einallv. lshen et alj ( l2014h and Ivan de Voort et'ffi] ^2015 ) esti¬ 
mated the enrichment history of r-process elements in the Galaxy, 
as traced by the [Eu/Fe] ratio, using a high resolution cosmologi¬ 
cal zoom-in simulation. Unlike previous studies, it was found that 
the nucleosynthetic products from compact binary mergers can be 
incorporated into stars of very low metallicity and at early times, 
even with a minimum time delay of 100 Myr and that compact 
binary mergers could be the dominant source of r-process nucle¬ 
osynthesis. 

In the present Section, we bring our own contribution to this 
difficult debate. We first summarize the available Eu observations 
and then propose a detailed comparison with our predictions of the 
Eu evolution along the cosmic history, using the three possible cos¬ 
mic SER density discussed in § 12.11 


4.1 Eu: observations 


We have gathered most of the recent observations on Eu abun¬ 
dance (Eu/H and [Eu/Ee] as a function of [Ee/H]) in metal 


poor 

stars jFrancois et a 

lJl2007l; 

Honda et alj|2004l: iBarklem et al. 

2005 

; Simmerer et al. 

12004: 

RoedereJ 201 ll; Roederer et al. 

201C 

. l20li 2014Jbllct 

iRenetalJ 20121: IWorlev et al. 20131 to- 


gether with the abundances found in dwarf spheroid a l sys- 


terns JShertone. Cote & Stetson 1 200 ll: Shertone. et al. 1 

2003; 

Geilser et al. 1 20051; Cohen and Huang 20091; Letarte et al. 1 

201C; 

Starkenburg et al. 112013; Me William et al. 1 2013). It is interesting 


to note that the bulk of data coming from external dwarf galaxies 
is embedded in galactic ones. The solar abundances are taken from 
iLodders I (l2003l) . 


The determination of Eu abundance suffers from the fact that 
only few transitions are visible in the spectrum of metal poor stars. 
The strongest transitions are rather weak and it becomes particu¬ 
larly difficult to detect and estimate the abundance of this element 
at low metallicities. Typically, the strongest Eu line becomes al¬ 
most undetectable below [Fe/H] -3.5. The detection or the non¬ 
detection depends on the S/N ratio of the spectrum, on the presence 
or not of an overabundance of Eu in a given star and finally on the 
temperature and gravity of the star. This situation clearly shows up 
in the diagrams where the number of true detection of Eu at very 
low metallicity are only a handful. Moreover, due to the scarcity of 
very metal poor stars in the halo, the new very metal poor candi¬ 
dates are faint and require a significant fraction of observing time 
on 10 m class telescopes. Consequently, the number of stars with 
[Fe/H] -3.5 for which the Eu abundance has been measured is 
small. The sample used in the present paper does not represent an 
unbiased sample of stars, so the distribution function of [Eu/Fe] vs 
[Fe/H] has no significance. In several stars, the upper limit of the 
Eu abundance has been computed. These upper limits are useful 
as they complement the real measurement and show that, so far, 
no star below [Fe/H] =: -3.5 has a very high [Eu/Fe] ratio, i.e. of 
the order of what is found in stars at [Fe/H] -3.5. Although we 
cannot fully rule out the existence of such stars with [Eu/Fe] >1.2 
dex, the present sample seems to indicate that the upper envelope 
of the [Eu/Fe] stops to rise at [Fe/H] ^ -3.0 This hypothes is re¬ 
ceives support from the recent study of iHansen et al.1 bOldl) who 
analyzed 4 stars with [Fe/H] < -4.00 and did not report any mea¬ 
surement of the Eu abundance although the wavelength range of 
the spectra used in this study covers the Eu transitions wavelength. 

The inspection of observed abundances (see Eig.|7]l shows the 
apparent existence of two branches starting below [Ee/H] ^ -2.5, 
the first one with high values of [Eu/Ee] as high as -1-2.0, the sec¬ 
ond one with decreasing values down to [Eu/Fe] == -1.0 at [Fe/H] 
-3.2. As the sample is biased, it is difficult to be fully confident 
about the reality of this two branches. The plot could also well be 
interpreted as an increase of the dispersion as the metallicity de¬ 
creases. Indeed, semi analytic models of merger tree which trace 
the chemical evolution of matter in different regions of the Uni¬ 
verse (from over- to under-dense regi ons) show clearly a d ispersion 
growing as a function of the redshift JOvorkin et al.ll201^ . 

The origin of the stars with a high level of r-process ele- 
ments like Eu is still a matter of debate. A recent analysis by 
iRoederer et al.l ( I2014bl l3l has shown and confirmed that this over¬ 
abundance is also found in main-sequence turn-off stars rejecting 
the hypothesis that the star would have produced itself such over¬ 
abundances. They also did not find any compelling evidence to 
suggest that a noticeable high fraction of highly r-process enriched 
stars are members of binary systems, assigning the origin of the r- 
process enhancement to a companion star. They also confirmed that 
these peculiar stars do not present a different fr-element chemical 
signature from the bulk of the other metal poor stars. Therefore, the 
site responsible for the production of this r-process enhancement is 
not expected to produce any chemical anomalies for light elements. 
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Figure 7. Cosmic evolution of Eu: comparison between the two possible astrophysical sites Eu/H and [Eu/Fe] vs [Fe/H]. Evolution of Eu/H and [Eu/Fe] 
as a function of [Fe/H] either in the CCSN (blue lines) or in the NSM (black lines) scenario. The evolution is computed for the three SFR modes considered 
in the present study, SERI (solid lines), SER2 (dotted lines), SFR3 (dashed lines). In the CCSN scenario, the Eu yield is 10“^ Mq per supernova. In the NSM 
scenario, the Eu yield is 7 X 10“^ Mq per merger, the fraction of binary compact objects is a = 0.002 and the coales cence timescale is AInsm = 0.2 Gyr. Data 
points come from diffe rent reference s: black points and u pper limits at lo w metallicity fro m ^Francois et alJl2 007ll. brown points at higher metallicity from 
'simmereretalj20M ), yellow from jfiarklem et al.l200^ . magenta from iRen et alj2012h . red from jRoederer et alj201^|2014a|l 3), heavy blue points from 


Roederer t 


). The bulk of black points at high metallicity come from dwarf spheroidal systems jSherton^*C6td & Stetson ll200ltrShertone. et al. I 


[^OrotlGeilseretal. l2005l: ICohen and Huang l2009l: iLetaife et al. I20l3: IStarkenburg et al. l2013l: iMcWilliam et al. 1201^ 



[Fe/H] [Fe/H] 


Figure 8. Same as Fig.|2]but without error bars. 


4.2 Eu: yields from CCSN and NSM 

Each of the two possible r-process sites, CCSN and NSM, may po¬ 
tentially contribute to the galactic r-process enrichment, but we will 
assume in the present work that they are mutually exclusive, so that 
the the r-process in the Universe is either originating from CCSN or 
NSM but not from a mix of both sites. While the cosmic evolution 


of r-process elements in the CCSN scenario essentially depends on 
the amount of r-process material ejected by the explosion, in the 
NSM scenario, it is also function of the binary fraction of compact 
objects a and the coalescence timescale of the binary system AInsm- 
We describe here the adopted yields of Eu expected to be ejected 
from CCSN or NSM. 
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4.2.1 Eu yields from CCSN 


Since there exists so far no consistent r-process calculation in 
CCSN, the Eu yields from supemovae remain highly uncertain. The 
CCSN Eu yields are therefore free parameters that can only he ad¬ 
justed to reproduce some observables, like the [Eu/Ee] trend as a 
function of the [Fe/H] observed in galactic stars or the total amount 
of Eu material in the Galaxy. In particular in this latter case, if we 
assume that the heavy-element composition in the solar system is 
representative for the whole Galaxy (as suggested by the small scat¬ 
ter of the Eu abundance in t he present-day Milky Way), for a Eu 
mass fraction of 3.7 x 10“'° jLodders II2003I : lAsplundjetalJ 200^ 
and a total galactic baryonic mass of 6x lO*" ( lMcMillanl201 ih . 
this implies a Eu Galactic content of 22 Mq. Assuming that the 
CCSN events occur at a constant rate identical to the currently ob¬ 
served one, i.e. about one per century during the approximate age 
of the Galaxy of 1.3 x 10'° yr, some 1.3 10® type-II explosions 
should be responsible for the 22 Mq of Eu, i.e. each CCSN should 
in this case produce about 1.7 x 10“^ of Eu. For this reason, by 
simplicity, a constant yield of 10“^ Mq of Eu is adopted here for all 
CCSN, regardless of the mass of the progenitor. 


4.2.2 Eu yields from NSM 

While a successful r-process cannot be obtained ab initio in CCSN 
simulations, NSM models consistently predict a significant pro¬ 
duction and ejection of r-process material. In particular, a mass- 
weighted combination of the dynamical ejecta from both the bi¬ 
nary merger phase and the secular eject a from the BH-torus evo¬ 
lution has been shown jjust et al.ll2015h to reproduce the solar r- 
abundance pattern, and therefore also the elemental one seen in 
ultra-metal-poor stars, amazingly well in the whole 90 £ A < 210 
range. This calculation represents today the most consistent es¬ 
timate of Eu production by NS-NS and NS-BH mergers taking 
consistently into account both the prompt ejecta as well as the 
relic BH-torus outflow contributions. For this reason, the Eu yields 
obtained in such combined models are considered in the present 
study. More precisely, the three combined systems (merger model 
plus remnant model) corresponding to model s with torus ma sses 
Aftorus = 0.03,0.1 and 0.3 (see Fig. 19 of ijust et alj|2015h are 
found to eject between 7 x 10“^ and 2 x lO^'* Mq of Eu. A fiducial 
Eu yield of 7 x 10“^ Mq is adopted here ( see the distribution of the 
ejected mass in NSM mergers in Fig. 6 in iFrver et al. I ( 1201 Sh ). 


Regarding the NSM scenario, with the fixed value of 7 x 
10“^ Mq for the Eu yield (see § 14.2.2b . observations at [Fe/H] > 
-2.5 are seen to be in good agreement with the predicted evolution 
for a fraction of binary compact objects a = 0.002 and a coales¬ 
cence timescale AtNSM = 0.2 Gyr. A larger Eu yield of 2 x 10“'' Mq 
would logically imply a lower value of a =: 0.0007, both parameters 
being degenerate in this approach. The agreement with the observa¬ 
tions of Eu at low metallicity depends strongly on the coalescence 
timescale, which can therefore be constrained for each SFR mode, 
as illustrated in Figs. muni and where we compare four coa¬ 
lescence time scales AtNSM = 0, 0.05, 0.1, 0.2 Gyr for each SFR 
mode. The best fit is obtained for a typical coalescence timescale 
of 0.2 Gyr for SERI (Fig.|9j, 0.1 Gyr for SFR2 (Fig.[T0ll and 0.05 
Gyr for SFR3 (Fig.[nj. 

The dilference stems from the [Fe/H] evolution which differs 
significantly for the three SFR modes, as already pointed out in 
§ 13.2.21 (see Fig. |5j. The more iron produced at early time, the 
shorter coalescence timescale needed to explain the low-metallicity 
r-process-rich stars. An important conclusion, illustrated by these 
results, is that a coalescence timescale larger than typically 0.3 Gyr 
is excluded by the observations of Eu at low [Fe/H] for all SFR. 
Note that the production of iron from type la supernovae is in¬ 
cluded in our calculation but their contribution becomes efficient 
only at the end of the cosmic evolution, i.e. for [Fe/H] > -1. In¬ 
deed, due to the lifetime of the 2 Mq stars (2 Gyr) and the potential 
time delay of explosions, SNIa iron enriches the ISM at z < 3. It 
should also be noted that, contrary to the CCSN case, the Pop III 
component in the SFRl mode does not affect the predictions in the 
NSM scenario due to the time delay of the merging. 


4.3.2 Evolution of Eu as a function of redshift 

Fig.ll2lshows the Eu evolution as a function of redshift for the three 
SFR modes. To include observational data as a function of z, each 
[Fe/H] value needs to be associated with a corresponding redshift. 
In our approach, the iron abundance increases monotonically with 
time, so that there exists a unique relation between [Fe/H] and z, us¬ 
ing equation 1, but this relation is however function of the adopted 
SFR since the Fe evolution depends on the SFR, as shown in Fig.|5] 
(lower panel). For this reason, the three modes are shown in differ¬ 
ent panels. In the NSM scenario, the coalescence timescale adopted 
for each SFR mode is AfjvsM =0.2,0.1,0.05 Gyr for SFRl, SFR2, 
SFR3, respectively, as discussed above. 


4.3 Eu: cosmic evolution 

4.3.1 Evolution of Eu as a function of [Fe/H] 

The Eu cosmic evolution is computed for the three SFR modes. 
Figs. HI and [8] (with and without error bars) show the predicted 
evolution of Eu/H and [Eu/Fe] as a function of [Fe/H] for both 
r-process sites, separately, the CCSN (blue lines) and NSM (black 
lines), with a comparison to the observations described in Sect. l4.ll 
Concerning the CCSN scenario, the Eu yield of 10“’ Mq is 
seen to reproduce fairly well the spectroscopic data for [Fe/H] > 
-2.5. At lower metallicity, Eu is clearly overproduced by CCSN, 
specifically if we consider Eu data at lowest metallicity, as pointed 
in Sect. 14.11 In the SFRl mode, where a Pop III contribution is 
included, the overproduction is even more significant due to the 
lack of iron produc t ion at very low metallicities in the models of 
IWooslev & Weav^ l ll995h (Case B), whereas Eu is assumed to be 
produced in all mass ranges of massive stars. 


4.4 Sensitivity analysis 

4.4.1 Eufrom CCSN and the Pop III star component 

As already mentioned, the Eu nucleosynthesis by CCSN remains 
highly uncertain. In our standard calculations (§ 14.3t , it was as¬ 
sumed that each CCSN produces 10“’ Mq of Eu regardless of the 
mass and metallicity of its progenitor. If we consider in the SFRl 
mode the existence of a Pop III star component, as required to ex¬ 
plain the Thomson optical depth of the CMB measured by WMAP9 
(see Fig. 0, Fig. H] clearly shows a large overproduction of Eu at 
low metallicity. So, would the r-process successfully take place in 
CCSN, the observation of Eu in ultra-metal-poor star could be a 
discriminant indicator highlighting the existence or not of Pop III 
stars. However, if we assume that Pop III stars, or more generally 
Z < 10“''Zq stars, do not produce r-process elements significantly 
(less than typically 5 10“° Mq), the Eu evolution curve in SFRl 
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Figure 9. Cosmic evolution of Eu in the NSM scenario: effect of the coalescence timescale (SFRl case). Same as Fig.[7]in the NSM scenaiio and SFRl 
case. We compare four different coalescence timescales: 0 (dot-dash line), 0.05 (dotted line), 0.1 (dashed line) and 0.2 Gyr (solid line). 



H 



Figure 10. Cosmic evolution of Eu in the NSM scenario: effect of the coalescence timescale (SFR2 case). Same as Fig. |9]in SFR2 ( 



Figure 11. Cosmic evolution of Eu in the NSM scenario: effect of the coalescence timescale (SFR3 case). Same as Fig.l^in SFR3 case. 
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Figure 12. Cosmic evolution of Eu: comparison between the two possi¬ 
ble astrophysical sites (2) En/H vs redshift. Evolution of Eu/H as a func¬ 
tion of z either in the CCSN (blue lines) or in the NSM (red lines) sce¬ 
nario. The evolution is computed for the three SFR modes considered in 
the present study, SFRl (upper panel), SFR2 (middle), SFR3 (lower panel). 
The three SFR modes are shown in different panels because the observa¬ 
tions expressed in terms of [Fe/H] do not conespond to the same redshift z 
when different SFR are considered. The yields in the two scenarios are the 
same as in Fig.|2] In the NSM scenario, the binary fraction is cr = 0.002 and 
the coalescence timescale is 0.2 Gyr for the SFRl mode, and 0.1 Gyr for 
the SFR2 and 0.05 Gyr for SFR3. 


model, becomes almost equivalent to the predictions obtained with 
SFR2 and SFR3, as shown in Fig.fO] 


4.4.2 Sensitivity to the early iron production 

Our results depend on the adopted scenario for star formation in 
the Universe. This is why a detailed sensitivity study has been per¬ 
formed with different SFR stories, different stellar yields, with the 



[Fe/H] 


Figure 13. Cosmic evolution of Eu in the CCSN scenario: effect of the 
Eu yield (SFRl case). Evolution of Eu/H as a function of [Fe/H] in the 
CCSN scenai'io for the SFRl mode, assuming either the standard prescrip¬ 
tion for the Eu yield (solid line), i.e. 10“^ M© per supernova for massive 
stars at all metallicities, or a modified prescription (dashed line) where 
Eu is produced in supernovae only for massive stars with a metallicity 
Z > 10-^Z„. 


aim of fitting the maximum of observational constraints, i.e. cos¬ 
mological and local. _ 

Recent studies publi shed by other authors jMatteucci et al.l 
I 2 OI 4 I : iKomiva et ^l2014h concluded that a merger time scale of 
the order of 1 - 10 Myr was necessary to reproduce the observed 
evolution of r elements with NS mergers alone. We find a different 
conclusion constraining the merger time scale to the range of 100 
- 200 Myr. We made a careful analysis of the origin of this dis¬ 
crepancy. We show that the merger timescale depends strongly on 
the iron evolution which itself depends on the global astration rate 
and the choice of yields in supemovae or in hypernovae. Unfor¬ 
tunately, the CCSN production of Fe is poorly known. It depends 
essentially on the explosion energy and on the mass cut (fraction 
of the mass trapped in the remnant). Iron ejection is more im¬ 
portant for higher en ergies of the explosi o n. Typ ically, we choose 
the SNII models by IWooslev & Weav^ ( 119951) = 10^' 

ergs). At solar metallicity and for stars of 40 (20) Mq, they give 
an iron yield of 0.03 (0.02) Mq respectively and 0.02 (0.1) Mq 
res pectively at Z = 0.001. On the other hand, hypernova models 
by iKobavashi et aP ) l2006h {Eexpiod. = 10 - 30 x 10^ ‘ ergs) cor¬ 
respond to 0.277 (0.035) Mq of iron for the same stars at solar 
metallicity and 0.26 (0.08) Mgrespectively at Z = 0.001. Indeed, 
hypemovae are much more energetical than the SNae. As seen 
in Fig. M a higher Fe early production constraints the merge r 
time scale to be sho rter. The studies by iMatteucci et al.l (|2014|): 
iKomiya et al.l (1201 4l) assume stellar yields from Kobavashi et al.l 


( 200^ (hypernova case), w hereas we assume the stellar yields from 
IWooslev & Weaved 1995h (supernova case). We favor this assump¬ 
tion, since obser vations of the most meta l-poor stars show a low 
iron abundance dFrebel and Norrisll2015h . It is worth noting that 
this result leading to merger timescales of the order of 100 - 200 
Myr is in better agreement with our knowledge of the properties of 
the binary NS population (see Sect. l2.2.3t . 
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Figure 14. Cosmic evolution of Eu in the NSM scenario: sensitivity to 
the iron production and to the lower end of the IMF. Upper panel: evolu¬ 
tion of [Fe/H] as a function of time in the SFR2 case, c omputed either with 
our st andard prescription (solid line: Fe yields from IWooslev & Weav^ 
(HH), case B, and Mi„f = 0.1 Mq f or the IMF), or w i th a m odified pre¬ 
scription (dashed line: Fe yields from iKobavashi et alj fcOOdl . hypemova 
case, and M\^f = 0.8 Mq for the IMF) ; Lower panel: corresponding evo¬ 
lution of [Eu/Fe] as a function of [Fe/H]: our standard prescription for iron 
production (solid line) with a coalescence timescale of A(nsm =0.1 Gyr is 
compared with the modified prescription for three values of the time delays, 
i.e. 0.1 Gyr (dashed line), 0.01 Gyr (dotted line), 0.001 Gyr, (dot-dash line). 


100 Mq range, obtain the same results regarding the production 
of iron. However, as s aid above, other studie s find different con¬ 
clusions; i n particu lar, iKobayashi et ^ ( 120061) predict higher iron 
yields than IWooslev & Weaved i ToOSh for massive stars regardless 
of the metallicity. Fig. [14] illustrates the impact of Fe yields on 
the Eu evolution when considering the yields of iKohavashi et aU 
(l2006h . Note that we have also modified the lower end (Minf) of the 
IMF from 0.1 to 0.8 Mq to increase the gas processing in stars and 

favor again the iron pro duction. _ 

With the yields of IKohavashi et all ( l2006h . iron at the very 
beginning of the evolution is seen to be synthesized earlier than in 
our standard model, leading to an overall abundance more than 10 
times higher (Fig. [m upper panel). We compare in Fig.ll4l (lower 
panel) our standard SFR2 model with A?nsm = 0.1 Gyr (solid black 
line) w ith the results obtained with the yields of IKohavashi et aU 
( l2006h and three different coalescence timescales, namely 0.1 Gyr 
(dashed line), 0.01 Gyr (dotted line) and 0.001 Gyr (dot-dash line). 
Consequently, a significant reduction of the coalescence timescale 
(AfNSM = 1 Myr) needs to be applied to compensate for the early 
increase in the Fe production in this conte xt. It should be recalled 
that it is close to the scenario obtained by iMatte ucci et aP ( l2014h 
using the same Fe yields from IKohavashi et alj i 20061) and an IMF 
with Minf = 0.8 Mq. Clearly, the more iron is produced at early 
time, the shorter the predicted coalescence timescale required to 
explain the Eu abundance in the ultra-metal-poor stars. 


4.4.3 Time delay distribution and dispersion in data 

When fitting the observed early cosmic evolution of Eu in the 
NSM scenario, the coalescence timescale AtNSM is constrained to 
be short, typically in the range 50-200 Myr depending on the SER 
(see § 14.3. H i, as binary systems with longer timescales are not 
merging at large redshifts. As discussed in Sect. 12.2.31 such sys¬ 
tems with a relatively low value of Af^sM represent a significant 
fraction of all the systems, but not all of them. In Eig.[T5| we show 
the evolution of Eu using a broader, and more realistic, distribu¬ 
tion of coalescence timescales. Based on the distribution plotted 
in Eig. 13 we assume 13% of the NSM in the 0-10 Myr range 
(AfNSM = 5 Myr), 10% in the 10-100 Myr range (AInsm = 30 Myr), 
15% in the 100-500 Myr range (AfNSM = 300 Myr), and 62% with a 
coalescence timescale above 500 Myr (AInsm = 1 Gyr)- The calcu¬ 
lation is carried out for the SERI mode and we present the result for 
both a = 0.002 (solid line) and a - 0.004 (dotted line). When com¬ 
pared to Fig. 13 the early evolution is slightly modified due to the 
contribution of systems with even shorter coalescence timescales 
than in our reference case. As mergers are rare events with large 
Eu yields (compared to CCSN which are more frequent but pro¬ 
duce less Eu), a large dispersion can be expected at low metallicity 
when only a small number of mergers affected the local environ¬ 
ment. This could explain the observed dispersion shown in Eig. [7] 
and discussed in Sect. 14. II 


As mentioned, the cosmic Fe enrichment plays a fundamen¬ 
tal role in the interpretation of the chemical abundance evolu¬ 
tion. To analyse the impact of this input on the determination of 
the NSM coalescence timescale, the Fe yield produced by mas- 
si ve stars has been va r ied in the SFR2 case. The standard yields 
of IWooslev & Weaveil ( Il995l) . more specifically their case B (for 
which low-metallicity high-mass stars do not p roduce iron), have 
been used. Note that dHeger & WooslevI 120101) which follow the 
evolution and the explosion of zero-metallicity stars in the 10 - 


5 THE PREDICTED NSM RATE 

Once the fraction a and the coalescence timescale AInsm have been 
adjusted to reproduce the cosmic evolution of Eu, as described in 
the previous section, the cosmic evolution of the merger rate is 
determined as all free parameters entering in Eq. 0 are known. 
The result is plotted as a function of the redshift in Fig. |3 for the 
SFRl and SFR3 modes, together with the CCSN rate and the NS 
and BH birth rates. The SFR evolution is shown for comparison. 
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Table 1. Predictions for the merger rate within the horizon of advanced Virgo/LIGO obtained with the constraint from the cosmic evolution of Eu using SFR1, 
SFR2, SFR3 models. For each SFR, the lower and higher values of the predicted rate coiTespond to the higher and lo wer limits on the yield of Eu in mergers, 
i.e. 2 X 10 ""^ Mq and 7 x 10~^ M q, respectively. The range of independent predictions compiled bv lAbadie et alj j201Ql) is given on the first line for comparison. 
Following lAbadie et alJi2Ql3) . we adopt 200 Mpc (resp. 420 Mpc) for the size of the horizon for NS-NS (resp. NS-BH) mergers. 


NS-NS merger rate (yr NS-BH merger rate (yr 


Abadie et al. (2010) 

40 (0.4-A00) 

10 (0.2-300) 

SFRl 

2.4-6.7 

2.7-7.7 

SFR2 

2.-5.7 

2.3-6.87 

SFR3 

3.8-10.9 

4.3-12.4 


Table 2. Predict ions for the kilonova rate for three different surveys: Palomar Transient Factory (PTF, limiting magnitud e of 21 and area of 2700deg^, 
iLaw et al. l20ld) . Large Syn optic Survey Te lescope (LSST, limiting magnitude of 24.5, area of 18000 deg^ . lAbell et al J 2 QO 9 I) and Euclid (limiting magnitude 
of 24.5, ai'ea of 15000 deg^. lLaureiis et al. IlioTl for SFRl, SFR2, SFR3 models. For each SFR, we consider either the most pessimistic case (Tow’) or the 
most optimistic case (’high’) corresponding to the highest and lowest possibles values of the yield of Eu, as in Tabled For each case, the range of predicted 
rates corresponds to the lowest and highest values of the peak luminosity of a kilonova (see text). 


Kilonova rate (yr *) 



PTF 

LSST 

Euclid 

SFRl low 
SFRl high 

0.0018-0.034 

0.005 - 0.096 

1.4-22.9 

4.1-65.4 

1.2- 19.1 

3.3- 54.5 

SFR2 low 
SFR2 high 

0.0014-0.028 
0.004 - 0.08 

1.2-20.1 

3.5-57.4 

1.0-16.7 

2.9-47.9 

SFR3 low 
SFR3 high 

0.003 - 0.054 
0.008-0.16 

2.3 - 36.4 
6.6-103.9 

1.9-30.3 

5.5-86.6 



-4 -3 -2 -1 


[Fe/H] 


Figure 15. Eu evolution with a time delay distribution. Same as in Fig.|3 

(SFRl case) for a distribution of coalescence timescales corresponding to 
Fig.|3](see text for more details). Solid and dotted curves correspond to two 
NSM fractions, of 0.002 and 0.004, respectively. 


5.1 Merger rate in the horizon of gravitational wave 
detectors 

On the basis of the cosmic evolution of the merger rate, it is possi¬ 
ble to estimate the rate expected in the horizon of the gravitational 
wave detectors advanced Virgo and ad LIGO and to c ompare those 
rates w ith independent predictions, as compiled by lAbadie et al.l 
(l2010l) . The result is plotted in Fig.[T^for the SFRl scenario. When 
considering the production of Eu in the NSM scenario, both param¬ 
eters, the Eu yield and the fraction of binary mergers a, are degen¬ 
erate when adjusted on the cosmic evolution of Eu. This is not the 
case for the calculation of the event rates. Therefore, we present in 
Fig. [T^ the results for two sets of parameters: either our standard 
Eu yield of 7 x 10“^ Mq and a binary fraction a = 0.002 (solid 
lines), or a larger Eu yield of 2 x 10““* Mq (as found in NSM sim¬ 
ulations, see Sect. 14.2.2i and a lower binary fraction a = 0.0007 
(dashed lines), both leading to exactly the same cosmic evolution 
of Eu. The predicted rates are given in Table[T]for the two extreme 
sets of parameters. 

Interestingly, the integrated merger rate in the horizon of grav¬ 
itational wave detectors obtained in the present study using our 
standard value for the Eu yield from NSM (7 x 10~^ Mq) is close 
to the central value of the predictions compiled by lAbadie et al.l 
( I 2 OIOI) . In contrast the upper value of the Eu yield tends to disfavo r 
the upper half of the predicted range given bv lAbadie et al.l ( l2010l) . 
as clearly illustrated in Eig. [m 


We also plot the lo cal merger rate taken from the compilation by 
lAbadie et alj ( |201(]|) . 


5.2 Kilonova rate 

As discussed in the introduction, NS-NS/BH mergers can po¬ 
tentially have two main types of electromagnetic counterparts, 
short GRBs and kilonovae. Short GRBs are associated to the pro- 
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Redshift 

0.05 0.1 0.15 




15 20 25 

Limiting magnitude 


Figure 16. The integrated merger rate (NS-NS, red line, NS-BH, blue 
line, and total, black line) as a function of the distance (lower axis) and 
the redshift (upper axis). Two vertical lines indicate the size of the horizon 
of the gravitational wave detectors advanced Virgo/LIGO. For comparison, 
the centr al value and range o f the predicted rates within this horizon com¬ 
piled by lAbadie et alj ilQlOl) is also indicated for both types of mergers. 
The calculation corresponds to the SFRl mode with the standard prescrip¬ 
tion (solid line: Eu yield of 7 x 10“^ M© and binary fraction a = 0.002) 
and with a modified prescription (dashed line: 2 x 10“"^ M© and binary frac¬ 
tion a = 0.0007), both leading to exactly the same cosmic evolution of Eu 
plotted in Fig. [7] 


duction of an ultra-relativistic outflow along the polar axis of 
the post-merger BH-torus system, with large uncertainties on the 
beaming angle of the jet. Therefore, the prediction of the short 
GR B rate from the merger rate requires a specific treatment (see 
e.g.lGuetta & Piran]l200d: IWanderman & Piranll20151: iFrver et aTI 

mm which we do not include in the present study. On the 
other hand, kilonovae are associated with a quasi-isotropic non- 
relativistic ejection that should be present in all systems. The lu¬ 
minosity of the kilonovae is powered by the radioactive decay of 
freshly synthesized elements and is consequently directly related to 
the nucleosynthesis of r-process nuclei. Note that NS/NS mergers 
are better sources of electromagnetic emission, as the ejected mass 
is very low in a noticeabl e fraction of NS/BH mergers (see Fig. 7 in 
iBelczvnski et alj l l2008bh l. As NS/BH mergers are about one order 
of magnitude less frequent than NS/NS mergers (see, Fig. [Tbt . this 
does not affect our results. 

Therefore, the detection of a kilonova would represent a major 
step to probe the contribution of NS-NS/BH mergers to the synthe - 
sis of heavy elements in the Universe (see also lMetzger et al.l2015h . 
There is only one or two tentative detections to date. The first can¬ 
didateis in association w ith the short GRB 130603B dBerger et al. I 
l2013l : lThnvir et alj|2013h . However, this observation is not detailed 
enough (only a single photometric measurement) to prove firmly 
its kilonova nature. The second candidate is associated to GRB 
060614, whose short or long nature is am biguous, but is also based 
on a photometric near-infrared bump only jYang et al. l2015h . More 
detections, with a photometric and spectroscopic follow-up, will be 
required to progress further. The peak of the kilonova is expected in 


Figure 17. The integrated kilonova rate per deg^ as a function of the lim¬ 
iting magnitude for the two same cases as in Fig. i.e. SFRl with the 
standai'd prescription (blue) or the modified prescription (red) for the Eu 
yield. The width of the shaded area in each case takes into account the the¬ 
oretical uncertainties on the peak luminosity of the kilonova (see Eq. jl2H . 
For comparison, horizontal yellow and green lines indicate the limit of 1 
event per year for three surveys (PTF, LSST and Euclid, see Sect. 15.21 and 
Table|2]for the limiting magnitude and area of each survey). 


the near-infrared range a few d ays after the m erger (it may be bluer 
in the case of NS/BH mergers, fa et al. Il2014h . The long-term evo¬ 
lution of the ejecta, which is decelerated by the ambient medium, 
also leads to radio emission over weeks. This signal is not con¬ 
sidered here though it may also become detectab le by future radio 
surveys (iNakar & Piran 1201 lllHran et al. l2013h . 

Our model allows to constrain the expected rate of kilono¬ 
vae for a given large survey as a function of the limit magnitud e 
in the optical/near-infrared range. Following lMetzger et alj ( l201(lf) . 
the peak luminosity scales as 


i'KN.peak - 5 X 10'“’ erg/s 


/ W V 

10-6 0.1c/ \10“^^o 


1/2 


( 12 ) 


where v and Mej are the velocit y and mass of the ejecta respectively 
and f is the heating efficiency JAmettll 982l:lLi & Paczvn skillQQsl : 
iMetzger et al.ll20ig) . Compared to Metzger et al. | |201(]|) . the nor- 
malization in Eq.[T2]has been corrected for a more realistic opac¬ 
ity of the r-process elements, which is expected to be ~ 10^ times 
larger compared to the opacity of iro n group e lements considere d 
in previous studies jKasen et alJ2013h (see also lKasen et alJ201.'5h . 
Using the extreme values obtained in the systematic s tudy o f the 
dynamics of the ejecta carried out bv iBauswein et al 1 20131). this 
leads to LxN.peak - (2.0 - 14) x 10'*° erg/s tsee also lGrossman et al. I 
|2014 . 


Assuming that all kilonovae have similar peak luminosities 
LxN.peak. we Can deduce from the merger rate the number of kilo¬ 
novae per year above a given flux. This leads to the results plotted 
in Fig. [T3 where the rate of kilonovae per year and per square de¬ 
gree is plotted as a function of the limiting magnitude of a given 
survey, in the SFRl case with the two same sets of parameters 
(Eu yield, compact binary fraction a) as in Eig. [761 The width of 
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the shaded area takes into account the range of ixN.peak discussed 
above. Horizontal lines indicate the limit above which at least one 
event per year is expected for different surveys, taking into account 
their magnitude limit and their field of view. The results are listed in 
Tablel^for each SFR and for the two extreme values of the yield of 
Eu. It appears that the current Palomar Transient Factory (PTE) sur¬ 
vey (limiting magnitude of 21 and area of 2700 deg^) fails to detect 
kilonovae, whereas the predictions are more optimistic for future 
surveys like the Large Synoptic Survey Telescope (LSST, limiting 
magnitude of 24.5, area of 18000 deg^) and mostly Euclid (limiting 
magnitude of 24.5, area of view of 15000 deg^) with a well adapted 
spectral range, for which typically ~ 2-100 events per year can be 
expected. 


6 CONCLUSIONS 

We have computed the predicted evolution of the cosmic rates of 
CCSN and NSM for three possible history of the SFR density, 
which are representative of the current uncertainties at high red- 
shift and reproduce important observational constraints such as the 
cosmic evolution of iron or the Thomson optical depth to CMB. We 
have deduced the corresponding expected evolution of Eu, a typi¬ 
cal r-process elements, considering that the astrophysical sites for 
the r-process are either CCSN or NSM. In the latter case, the com¬ 
parison to observations is used to constraint two free parameters, 
the binary fraction of compact objects and the typical coalescence 
timescale in NS binary systems. We obtain the following results: 

(i) The r-process evolution using the NSM scenario as the main 
astrophysical site is in good agreement with observations, assuming 
that the early evolution is dominated by mergers of binary systems 
with a coalescence timescale of the order of AJnsm ~ 100 Myr. 
Such mergers represent a significant fraction of all mergers accord¬ 
ing to recent estimations obtained with detailed population synthe¬ 
sis codes; 

(ii) In this merger scenario, the explosive events responsible for 
the production of r-process elements are more rare than in the 
supernova scenario, but each event yields more Eu. Therefore, a 
larger scattering in the r-process abundance is expected in metal- 
poor stars enriched by mergers. We present calculations of the early 
evolution of Eu obtained when summing up contributions of merg¬ 
ers with a broad range of coalescence timescales according to pre¬ 
dictions of population synthesis models and we find a reasonable 
agreement with the level of scatter observed. At high metallicities, 
the cumulative effect leads to a low sensitivity with respect to the 
coalescence timescale; 

(iii) The precise constraint on the dominant coalescence timescale 
relevant for the early evolution of Eu depends on the adopted SFR 
at high redshift, for which we have tested three scenarios repre¬ 
sentative of the current uncertainties, and on the iron yields of 
stars at low metallicity, for which we have compared two th eo¬ 
retical predictions ma de by either IWooslev & Weav^ ( Il995h or 
iKobavashi et al.l i2006h . The conclusion is that shorter coalescence 
timescales are favored by scenarios for which the production of iron 
at early times is more intense, i.e. for larger SFR at high redshift and 
larger iron yields coming from CCSN. In the standard conditions, 
we constrain the coalescence timescale Ar^sM to be in the range of 
50-200 Myr; 

(iv) In the merger scenario which is favored by observations at low 
metallicity, two additional parameters are degenerate: the fraction 
O' of NS in a binary system with another NS or a BH and the yield 
of Eu per merger. This second parameter is constrained to be in 


the range 7 x 10“^-2 x 10^“* Mq per merger in the most recent nu¬ 
cleosynthesis calculations. Adopting these two values leads to an 
upper and a lower limit on a (0.002 - 0.0007) which constrain the 
merger rate as a function of the redshift in the Universe; 

(v) We use the previous calculation to predict the expected rate 
of mergers within the horizon of gravitational wave detectors. We 
find that the predicted rate as constrained by the cosmic evolution of 
Eu falls well within th e range of prediction s by other independent 
methods compiled by lAbadie et alj l l201(]tl but clearly favors the 
mid or lower values, and not the most optimistic predictions. We 
predict typically 2 to 10 NS-NS mergers per year and an equivalent 
number of NS-BH mergers per year; 

(vi) Considering the most recent estimates of the opacity of r- 
process elements, we also predicted the rate of kilonovae in cur¬ 
rent and future optical/near-infrared large surveys. Our results are 
compatible with the lack of kilonova detections by PTF and are 
optimistic for the future LSST and Euclid projects, though with a 
large uncertainty; from 1 kilonova per year in the most pessimistic 
case to 100 events per year in the most optimistic case. 

Our results can be improved in the future in several ways. 
Successful r-process in CCSN with well-defined Eu yields, includ¬ 
ing a possible dependence with metallicity, would be very useful 
to definitively check the consistency of this scenario with obser¬ 
vations. A better determination of the SFR at large redshift and a 
convergence of theoretical predictions for iron yields in primordial 
CCSN would also strongly reduce the uncertainties on the parame¬ 
ters of the model in the merger scenario. In addition, more measure¬ 
ments of the abundance of r-process elements in stars at even lower 
metallicity ([Fe/H] < -3.5) would also allow to better constrain the 
model parameters. Finally, the predictions for the kilonova detec¬ 
tion rate in future surveys could be improved by taking into account 
the correction for the spectral range of the instruments. Additional 
predictions for the rates of other transients associated with mergers, 
i.e. short gamma-ray bursts and radio flares, could also be deduced 
from the predicted merger rate to provide new diagnostics in the 
future. 
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